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Abstract 


We introduce a symplectic no-core configuration interaction (SpNCCI) frame- 
work for ab initio nuclear structure calculations, in a correlated many-body 
basis which encodes an approximate Sp(3,R) symmetry of the nucleus. Such 
a scheme potentially provides a means of restricting the many-body space to 
include only those highly-excited configurations which dominantly contribute 
to the nuclear wave function. We examine the symplectic symmetry structure 
arising in an illustrative ab initio SpNCCI calculation for °Li. We observe both 
the dominance of symplectic symmetry in individual wave functions and the 
emergence of families of states related by symplectic symmetry. 
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1 Introduction 


A long-standing goal of nuclear physics is to quantitatively predict the structure 
of nuclei and understand their excitation modes ab initio, i.e., directly from re- 
alistic internucleon interactions. However, in a traditional oscillator-basis no-core 
configuration interaction (NCCI) calculation [1,2], the dimension of the many-body 
basis explodes as the number of nucleons and included single-particle excitations 
is increased. The basis size which would be required in order to obtain quantita- 
tively accurate predictions for nuclei with more than just a few nucleons becomes 
prohibitively large [3, 4]. 

The symplectic no-core configuration interaction (SpNCCI) framework intro- 
duces a correlated many-body basis, one which encodes an approximate Sp(3, R) 
symmetry of the nucleus. Our aims in pursuing symplectic many-body symmetry 
are twofold: (1) to use symmetries to accelerate convergence of ab initio results 
and (2) to understand the symmetries underlying many-body correlations in nuclei, 
including emergent rotation [5]. 

The symplectic group in three dimensions Sp(3,R) [6] enters into the nuclear 
many-body problem both through its relation to kinematics (coordinates and mo- 
menta) and through its connection to the harmonic oscillator. It is generated by 
the bilinears in coordinates and momenta, i.e., operators of the form 2;2j;, rp;, 
pixj, and pp; (i,j = 1,2,3). It thus arises naturally in problems involving the 
coordinates and momenta. Symplectic symmetry has a close relation to collective 
deformation and rotational degrees of freedom [7], through its Elliott SU(3) [8,9] 
(and rotor model [10]) subgroup, which is generated by the orbital angular momen- 
tum and quadrupole operators. The group Sp(3, R) is also the dynamical group for 
the harmonic oscillator, which defines the underlying basis for nuclear configuration 
interaction (or interacting shell model) calculations. 

These connections permit the construction of the symplectic shell model formal- 
ism, originally proposed by Rosensteel and Rowe [11]. While the goal of “diago- 
nalising a realistic many-nucleon Hamiltonian in a Sp(3,R) D SU(3) shell model 
basis, to obtain a fully microscopic description of collective states from first princi- 
ples” was envisioned early on [12], only recently has sufficient progress been made 
in the computational framework for ab initio calculations to allow this goal to be 
revisited [13-16]. 

From the viewpoint of accelerating convergence in NCCI calculations, it is per- 
haps most important to note that the need for including highly excited configurations 
in NCCI calculations exists, in large part, because the kinetic energy induces strong 
coupling across shells. The kinetic energy operator, as a generator of Sp(3,R), con- 
serves Sp(3,R) symmetry. Thus, by reorganizing the many-body space according 
to irreducible representations (irreps) of Sp(3,R), it is broken into subspaces which 
cannot be connected by the kinetic part of the Hamiltonian, as portrayed in Fig. 1. 
Combining symplectic symmetry with the no-core configuration interaction frame- 


Symplectic No-Core Configuration Interaction Framework 19 


a6 a" Fa a 
an = 
ad mt ee 


Figure 1: Reorganization of the NCCI many-body space into symplectic irreps: 
(left) the kinetic energy connects many-body states related by 0 or +2 oscillator 
quanta (as indicated by arrows), but (right) only connects states within a symplectic 
irrep (shaded cones). The detailed structure of an Sp(3,R) irrep in terms of U(3) 
irreps N(A, ) is illustrated in the inset (bottom right). 


work potentially provides a means of identifying and restricting the basis to include 
only those highly excited configurations which dominantly contribute to the nuclear 
wavefunction, thereby (it is hoped) reducing the size of the basis necessary to obtain 
accurate results. 

In this contribution, we briefly present a framework for ab initio calculations of 
nuclear structure in a SpNCCI scheme (Sec. 2). We then examine the symmetries 
arising in a SpNCCI calculation for “Li (Sec. 3). We observe both the dominance of 
symplectic symmetry in the individual wave functions and the emergence of families 
of states related by symplectic symmetry. These results are from work described in 
Ref. [17]. 


2 Symplectic NCCI framework 


The SpNCCI basis for the nuclear many-body problem consists of states which 
reduce an Sp(3,R) D U(3) D> SO(3) subgroup chain, i.e., which are arranged into 
nested irreps of all three of these groups. A symplectic irrep is built starting from 
some U(3) irrep with some lowest number of oscillator quanta, the lowest grade irrep 
(LGI). Laddering with the symplectic raising operator (A « b'b'), which creates 
pairs of oscillator quanta, generates the rest of the symplectic irrep, which is an 
infinite tower of U(3) irreps of increasing oscillator number. In practice, for SoNCCI 
calculations, this irrep must be truncated at some finite number of oscillator quanta. 

To understand the quantum numbers for the SpNCCTI basis, we first note that 
the group U(3) = U(1) x SU(3) is the product of the U(1) generated by the harmonic 
oscillator number operator, with quantum number JN, and the Elliott SU(3) group, 
which provides quantum numbers (A, 11). Thus, a U(3) irrep is labeled by quantum 
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numbers w = N,(Aw, Hw). The structure of a symplectic irrep is uniquely defined 
by the U(3) quantum numbers o = No(Ac, Uo) of its LGI, which thus serve as the 
Sp(3,R) quantum numbers of the symplectic irrep. The detailed structure of an 
Sp(3, R) irrep in terms of U(3) irreps is illustrated in the inset of Fig. 1, specifically 
for the 16(2, 1) irrep of *Be.! 

At the bottom of the group chain, SO(3) is the orbital angular momentum group, 
which provides the quantum number L. While the Sp(3, R) subgroup chain describes 
orbital (spatial) structure, there is also the complementary SU(2) spin group, with 
quantum number S. Orbital and spin angular momenta couple to give total angular 
momentum J. The SpNCCI basis is therefore classified according to symmetry labels 
for the subgroup chain 


[Sp(3,R) > U(3) D $O(3)] x SUs(2) > SU,(2). (2.1) 
Ne Owe NyGuey Eh Ss J 


Recall that here N.,(Aw, Uw) are the U(3) quantum number of the basis state itself — 
so, if we are only interested in U(3) symmetry, we can simply designate these N (A, /1). 
Then, No(Aq, fic) are the U(3) quantum numbers of the LGI from which the entire 
symplectic irrep is built. 

When defining an NCCI basis, it is more transparent and convenient to work 
with the number of oscillator excitations Nex relative to the lowest Pauli-allowed os- 
cillator configuration for the nucleus. Configurations with N.x = 0 are traditionally 
known as “Ofw” configurations, those with Nex = 2 as “2hw” configurations, etc. 
Similarly, for a SpNCCI basis, we shall refer to the number N.,,¢x (or simply Nex) of 
excitation quanta for U(3) irreps and Nox for the LGIs of Sp(3, R) irreps. Thus, in 
the symplectic decomposition of the many-body space shown schematically Fig. 1 
(right), each disk (light shading) respresents states with a given Nex = 0, 2, or 4, 
while each cone (dark shading) represents a symplectic irrep with LGI at No ex = 0, 
2, or 4. 

In our SpNCCI framework, matrix elements of the Hamiltonian (or other opera- 
tors) in the symplectic many-body basis are determined from certain “seed” matrix 
elements by use of symplectic laddering operations. The approach is based on ideas 
initially proposed by Suzuki and Hecht [18,19], but these have now been extended 
to accommodate general, realistic internucleon interactions. 

In order to evaluate the matrix elements of an operator in the SpNCCI basis, we 
first decompose the operator in terms of SU(3)-coupled components, or unit tensors. 
We then apply commutation relations of these unit tensors with the symplectic 
laddering operators to recursively express all matrix elements in terms of matrix 
elements among certain lowest states, specifically, the LGIs. We explicitly expand 


'The U(1) label N appearing in the U(3) labels N(A, jz) is actually not the number of oscillator 
quanta, per se, but rather the dimensionless oscillator Hamiltonian, which also includes a zero-point 
contribution of 3/2 for each nucleon [12]. Thus, in the lowest Pauli-allowed oscillator configuration 
of *Be, the four s-shell nucleons contribute 3/2 each, and the four p-shell nucleons contribute 5/2 
each, giving N = 16, as in the 16(2,1) irrep shown. 
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Figure 2: Energy eigenvalues from an SpNCCI calculation for °Li. Energies are plot- 
ted us. J(J +1), as conventional for rotational analysis [5]. States below the dashed 
line are predominantly Ofw, while those above the line have predominantly 2hw 
(or higher) contributions. States having dominant U(3) x SU(2) contributions with 
Nex(A, #)S quantum numbers 0(2,0)1 and 2(4,0)1 are highlighted (solid symbols, 
as labeled). 


these LGIs in terms of the SU(3) no-core shell model [SU(3)-NCSM] basis and use 
the SU(3)-NCSM code LSU3shel11 [20] to evaluate the seed matrix elements. The 
recurrence then involves coefficients obtained through vector coherent state (VCS) 
methods [21-23]. 


3 Illustration: Symplectic symmetry in °Li 


Since SpNCCI calculations are carried out in an Sp(3,R) D U(3) basis, we can 
immediately extract the decomposition of each eigenfunction into contributions with 
different symmetry characters, for all the groups in the chain (2.1). That is, we 
can decompose the wave function according to any combination of the quantum 
numbers 0 = No(Ag, Uc) for Sp(3,R) and w = Ny(Aw, Mw) for U(3), as well as the 
more familiar orbital angular momentum LF and spin S. 

By examining such decompositions, we can start to answer some of the questions 
we laid out above: Do the highly-excited oscillator-basis contributions to the nuclear 
wave functions primarily come from low-lying symplectic irreps, suggesting that 
symplectic truncation may be feasible? And can an Sp(3,R) D U(3) irrep structure 
provide a useful classification scheme for understanding the nuclear eigenspectrum? 
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Figure 3: Decompositions of calculated eigenstates by oscillator excitation quanta 
Nex (left) and by the No.ex of the contributing symplectic irrep (right). Decomposi- 
tions are shown for the highlighted states from Fig. 2, with dominant contributions 
from U(3) x SU(2) symmetry 0(2,0)1 (bottom) or 2(4,0)1 (top). Only oscillator 
configurations with even values of Nex contribute to natural parity states. 


Here we consider a SpNCCI calculation of ®Li, for which the low-lying energy 
spectrum is shown in Fig. 2. This calculation is based on the JISP16 internucleon 
interaction [24], with no Coulomb interaction (7.e., pure isoscalar). For initial exam- 
ination and benchmarking purposes, we simply choose a space which is truncated 
at 6 oscillator quanta, with oscillator basis parameter hw = 20MeV. That is, we 
take all symplectic irreps with LGIs having up to Nomax = 6 excitation quanta, 
then truncate each symplectic irrep to basis states with up to Nmax = 6 excitation 
quanta. 

Each of the states in Fig. 2 is found to have at most one or two dominant 
U(3) x SU(2) contributions, that is, well-defined w and spin quantum numbers. 


*The resulting spectrum is identical to that for a traditional Nmax = 6 M-scheme NCCI calcu- 
lation, since both spaces include all intrinsic excitations of up to 6 quanta (this equivalence permits 
rigorous numerical validatation against traditional NCCI codes such as MFDn [25, 26]). While the 
traditional 1/-scheme NCCI basis consists of laboratory-frame Slater determinants (which include 
center-of-mass excitations and which, while having definite M/, are in general admixtures of all 
J > |M|), the SpNCCI basis is defined in the intrinsic frame (consisting of states which are free of 
center-of-mass excitations and which are also naturally J-coupled). Compare the resulting dimen- 
sion of 197 822 for the laboratory-frame M = 0 space [3] with 3484 for the intrinsic-frame J = 0 
space [27], for SLi at Nmax = 6. 
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However, we focus on two specific U(3) x SU(2) “families” of states in this discussion, 
identified by the filled symbols in Fig. 2. We shall see that these together form part 
of a larger family of states with the same dominant symplectic symmetry. 

Let us start by classifying the states in the spectrum of Fig. 2 according to 
contributions from configurations with different numbers of oscillator quanta. The 
calculated states below an excitation energy of ~ 20MeV, as demarcated by the 
horizontal dashed line in Fig. 2, are dominated by N.x = 0 contributions. That is, 
in traditional shell model parlance, they have predominantly Ohw character. Above 
this point in the spectrum, eigenstates dominated by Nex = 2 contributions, or 2hw 
states, begin to appear. 

To provide a more quantitative picture of this situation, decompositions by Nex 
are shown in Fig. 3 (left): specifically, decompositions for the Ohw states highlighted 
at the bottom of Fig. 2 are overlaid in the lower panel [Fig. 3(b)], and decompositions 
for the 2hw states highlighted at the top of Fig. 2 are overlaid in the upper panel 
[Fig. 3(a)].? Clearly the “Ofw” states are not described entirely in the valence shell, 
but rather are heavily “dressed” with excited oscillator configurations, nor are the 
“2hw” states entirely orthogonal to the valence space. However, the Ohw states are 
dominated by Nex = 0 contributions at the 60%-80% level, while, for the 2hw states, 
the N.x = 0 contributions fall below the 15% level. 

We now look at the full decomposition with respect to U(3) x SU(2) character, 
labeled by the U(3) x SU(2) quantum numbers Nex(A, 4), in Fig. 4. The states in 
the Ohw family highlighted at the bottom of Fig. 2 all have dominant U(3) x SU(2) 
contribution 0(2,0)1. [For the ground state of °Li, the U(3) structure — and 0(2,0)1 
dominance — was explored by Dytrych et al. [28], on the basis of ab initio SU(3)- 
NCSM calculations.}] Contributions from the other Ohw U(3) x SU(2) irreps appear 
at the few percent level. This pattern of U(3) x SU(2) contributions is shown for the 
1} ground state and 3} first excited state in Fig. 4(c,d). The states in the 2hw family 
highlighted at the top of Fig. 2 instead have dominant U(3) x SU(2) contribution 
2(4,0)1, as shown for the 1* and 3* states of that family in Fig. 4(a,b). 

While the U(3) x SU(2) decomposition begins to tell us about the symmetry 
properties of these states, it provides only partial information on the symplectic 
structure. Any SpNCCI basis state with Nox excitation quanta comes from a sym- 
plectic irrep with Noex < Nex. That is, the basis state may itself be part of an 
LGI (Nox = Nex), or it may be obtained by the action of the symplectic raising 
operator from a lower LGI (Nox < Nex). 

Each bar in the U(3) x SU(2) decomposition histograms in Fig. 4 represents 
the total probability contribution of many basis states from several U(3) x SU(2) 


3When interpreting these, and subsequent, decompositions in an oscillator basis, it is important 
to keep in mind that they reflect a particular choice of oscillator length (here, corresponding to 
hw = 20 MeV) and are subject to change as the oscillator basis is dilated to other values of hw. 
Therefore, one should be cautious in attaching undue physical significance, e.g., interpreting them 
to mean that any given state has fundamentally a “Ohw” or “2hw” nature with respect to some true 
mean-field vacuum. 
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Figure 4: Decompositions of calculated eigenstates by U(3) x SU(2) contributions 
Nex(A, 11).S: for the lowest 1* and 37 states (c-d) and for excited 1* and 3* states 
of predominantly 2hw character (a-b). The contributions are ordered, from left to 
right in the figure, first by Nex, and then by SU(3) and spin quantum numbers. 


irreps sharing the same quantum numbers. For instance, there are in fact 8 different 
U(3) x SU(2) irreps with quantum numbers 2(4,0)1, for °Li. One of these comes 
from the Nox = 0 symplectic irrep 0(2,0)1, while the rest are themselves LGIs of 
No,ex = 2 symplectic irreps. 

Before breaking each wave function down into contributions from specific Sp(3, R) 
quantum numbers, it is informative to simply look at the decomposition of each wave 
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Figure 5: Decompositions of calculated eigenstates by Sp(3, R) x SU(2) contributions 
Noex(Ac;Ua)S: for the lowest 1+ and 3* states (c-d) and for excited 1* and 3* 
states of predominantly 2hw character (ab). 


function with respect to how “excited” the contributing symplectic irreps are. De- 
compositions by No x are shown in Fig. 3 (right). 

Since the Nox = 0 contribution must be at least as large as the Nex = 0 con- 
tribution, clearly the Ohw states will have a significant Nox = 0 contribution. In 
particular, for the highlighted family of states with U(3) x SU(2) character 0(2,0)1 
[Fig. 3(d)], the No ex = 0 contributions are in the 70%-90% range. Thus, a sub- 
stantial portion of the excited oscillator contributions [Fig. 3(c)] actually comes 
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from Nox = 0 sympletic irreps. This is encouraging for the viability of symplectic 
truncation schemes. 

But what about the 2hw states? For the highlighted family of states with U(3) x 
SU(2) character 2(4,0)1 [Fig. 3(b)], recall that this 2(4,0)1 contribution could come 
either from the No x = 0 symplectic irrep 0(2,0)1 or from the 2(4,0)1 LGIs of 
Nox = 2 symplectic irreps. So, are these 2hw states dominated by contributions 
from Nox = 0 symplectic irreps or from Noex = 2 symplectic LGIs? From the No ex 
decompositions [Fig. 3(b)], it is immediately apparent that they are dominated by 
contributions from Nox = 0, at the > 70% level. 

The full decompositions with respect to Sp(3,R) x SU(2) character, labeled by 
the Sp(3,R) x SU(2) quantum numbers No ex(Ac, Mc)S, are shown in Fig. 5. The 
most notable feature apparent in Fig. 5 is the dominance of the contribution from 
the 0(2,0)1 symplectic irrep, for all these states. In fact, all the highlighted states in 
Fig. 2 — both the low-lying Ohw states with 0(2,0)1 U(3) x SU(2) character and the 
high-lying 2hw states with 2(4,0)1 U(3) x SU(2) character — receive their dominant 
contribution from this single Nox = 0 symplectic irrep, the 0(2,0)1 symplectic 
irrep. 

Thus, states with very different oscillator excitation content [Fig. 3 (left)], and 
consequently very different U(3) content (Fig. 4), nonetheless form a larger family 
of states sharing the same symplectic symmetry. This is fundamentally the idea 
of the classification scheme for the SpNCCI basis states in an Sp(3,R) D U(3) 
scheme (2.1): many U(3) irreps form a single Sp(3, R) irrep, as illustrated in Fig. 1 
(inset). However, here we find the same Sp(3,R) D U(3) organizational scheme 
holding for the physical spectrum of energy eigenstates obtained from an ab initio 
calculation. 


4 Conclusions 


In conclusion, a framework for ab initio no-core configuration interaction calculations 
in an Sp(3,R) D U(3) symplectic basis is now in place. This scheme builds on an 
existing SU(3)-coupled framework for the nuclear problem and combines it with the 
group-theoretical machinery for Sp(3, R). 

In initial calculations, illustrated here with °Li, we confirm Sp(3, IR) as an ap- 
proximate symmetry of states throughout the low-energy spectrum. These states 
are characterized by mixing of a few dominant symplectic irreps. Then, families of 
states arise with similar symplectic structure, despite their differing U(3) content. 

To take full advantage of approximate symplectic symmetry in nuclei, as a means 
of accelerating convergence of NCCI calculations, it will be necessary to go beyond 
the simple benchmark Na, truncation scheme for the SpNCCI space illustrated 
here. While some gains can likely be made by truncating the space to dominant 
No(Ao; Ma )S symplectic subspaces, these subspaces are highly degenerate, especially 
as we go to higher No ex. 
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The most effective truncation will therefore likely come by identifying the few 
dominant symplectic irreps from within these subspaces (e.g., by some variant of 
importance truncation [29]). Once these “Hamiltonian preferred” symplectic irreps 
are identified, a more stringent truncation can be carried out to a relatively small 
number of symplectic irreps within each subspace. 


Acknowledgements 


We thank David J. Rowe for invaluable assistance with the Sp(3,R) formalism and 
Chao Yang, Pieter Maris, Calvin W. Johnson, and Patrick J. Fasano for discussions 
of the computational implementation. 

This material is based upon work supported by the U.S. Department of Energy, 
Office of Science, Office of Nuclear Physics, under Award Number DE-FG02-95ER- 
40934, by the U.S. Department of Energy, Office of Science, Office of Workforce 
Development for Teachers and Scientists, Graduate Student Research (SCGSR) pro- 
gram, under Contract Number DE-AC05-060R23100, by the Research Corporation 
for Science Advancement, under a Cottrell Scholar Award, and by the Czech Science 
Foundation, under Grant Number 16-16772S. TRIUMF receives federal funding via 
a contribution agreement with the National Research Council of Canada. 

This research used computational resources of the University of Notre Dame Cen- 
ter for Research Computing and of the National Energy Research Scientific Comput- 
ing Center (NERSC), a U.S. Department of Energy, Office of Science, user facility 
supported under Contract DE-AC02-05CH11231. 


References 
1] P. Navratil, J. P. Vary, and B. R. Barrett, Phys. Rev. Lett. 84, 5728 (2000). 


2| B. R. Barrett, P. Navratil, and J. P. Vary, Prog. Part. Nucl. Phys. 69, 131 
(2013). 


3] P. Maris, J. P. Vary, and A. M. Shirokov, Phys. Rev. C 79, 014308 (2009). 

4| P. Maris and J. P. Vary, Int. J. Mod. Phys. E 22, 1330016 (2013). 

5] P. Maris, M. A. Caprio, and J. P. Vary, Phys. Rev. C 91, 014310 (2015). 

6] B. G. Wybourne, Classical Groups for Physicists (Wiley, New York, 1974). 

7| D. J. Rowe, A. E. McCoy, and M. A. Caprio, Physica Scripta 91, 033003 (2016). 


8] J. P. Elliott, Proc. R. Soc. London A 245, 128 (1958). 


9| M. Harvey, Adv. Nucl. Phys. 1, 67 (1968). 


28 


10 


11 
12 


13 


14 


15 


16 


17 


18 


19 


20 


[21] 


[22] 


[23] 


24 


25 


26 


27 


A. E. McCoy, M. A. Caprio, T. Dytrych 


J. Carvalho, R. Le Blanc, M. Vassanji, D. J. Rowe, and J. B. McGrory, Nucl. 
Phys. A 452, 240 (1986). 


G. Rosensteel and D. J. Rowe, Ann. Phys. (N.Y.) 126, 343 (1980). 
D. J. Rowe, Rep. Prog. Phys. 48, 1419 (1985). 


T. Dytrych, K. D. Sviratcheva, C. Bahri, J. P. Draayer, and J. P. Vary, Phys. 
Rey. Lett. 98, 162503 (2007). 


T. Dytrych, K. D. Sviratcheva, C. Bahri, J. P. Draayer, and J. P. Vary, Phys. 
Rev. C 76, 014315 (2007). 


T. Dytrych, K. D. Sviratcheva, C. Bahri, J. P. Draayer, and J. P. Vary, J. Phys. 
G 35, 095101 (2008). 


T. Dytrych, K. D. Sviratcheva, J. P. Draayer, C. Bahri, and J. P. Vary, J. Phys. 
G 35, 123101 (2008). 


A. E. McCoy, Ab initio multi-irrep symplectic no-core configuration interaction 
calculations, Ph.D. thesis, University of Notre Dame (2018). 


Y. Suzuki and K. T. Hecht, Nucl. Phys. A 455, 315 (1986). 
Y. Suzuki, Nucl. Phys. A 448, 395 (1986). 


T. Dytrych, P. Maris, K. D. Launey, J. P. Draayer, J. P. Vary, D. Langr, 
E. Saule, M. A. Caprio, U. Catalyurek, and M. Sosonkina, Comput. Phys. 
Commun. 207, 202 (2016). 


D. J. Rowe, G. Rosensteel, and R. Carr, J. Phys. A 17, L399 (1983). 
D. J. Rowe, J. Math. Phys. 25, 2662 (1984). 


K. T. Hecht, The vector coherent state method and its application to problems 
of higher symmetry, Lecture Notes in Physics Vol. 290 (Springer-Verlag, Berlin, 
1987). 


A. M. Shirokov, J. P. Vary, A. I. Mazur, and T. A. Weber, Phys. Lett. B 644, 
33 (2007). 


P. Maris, M. Sosonkina, J. P. Vary, E. Ng, and C. Yang, Procedia Comput. Sci. 
1, 97 (2010). 


H. M. Aktulga, C. Yang, E. G. Ng, P. Maris, and J. P. Vary, Concurrency 
Computat.: Pract. Exper. 26, 2631 (2013). 


F. Q. Luo, M. A. Caprio, and T. Dytrych, Nucl. Phys. A 897, 109 (2013). 


Symplectic No-Core Configuration Interaction Framework 29 


[28] T. Dytrych, K. D. Launey, J. P. Draayer, P. Maris, J. P. Vary, E. Saule, 
U. Catalyurek, M. Sosonkina, D. Langr, and M. A. Caprio, Phys. Rev. Lett. 
111, 252501 (2013). 


[29] R. Roth and P. Navratil, Phys. Rev. Lett. 99, 092501 (2007). 


